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^ : ABSTRACT 

o 

^j^, The deconfinement phase transition from hadronic matter to quark matter can con- 

^ I tinuously occur during spins down of neutron stars. It will lead to the release of latent 

I heat if the transition is the first-order one. We have investigated the energy release 

' of such deconfinement phase transition for rotating hybrid stars model which include 

mixed phase of hadronic matter and quark matter. The release of latent heat per baryon 

i-C ■ is calculated through studying a randomly process of infinitesimal compressing. Finally, 

I . we can self-consistently get the heating luminosity of deconfinement phase transition by 

O 

c/3 , stars. 



imputing the EOS of mixed phase, and based on the equation of rotation structure of 



O , 1. Introduction 

o 

^ I Neutron stars (NSs) provide us a unique playground to study the properties of super-dense 

' matter in the most extreme physical conditions. NSs participation in various astrophysical phenom- 

■ ena usually presented challenges for us. The different equation of state (EOS) and their properties 

■ of matter at high densities would play an important role in understanding some astrophysical 
phenomena. 



One of the most intriguing predictions of some theories of dense matter is a possibility of a 
phase transition into an 'exotic' state, including pion and kaon condensation, and deconfinement 
of quarks. Of them, the deconfinement phase transition represented the most profound effect on 
the structure and dynamics of NSs. Many investi gation s were interesting in th e phase transition 



which is of the first-order type(Pisalski & Wilczek (jl984l ) and Gavai et al. (|l987l )). In the simplest 
case, equilibrium phase transition from the normal, lower density phase to the pure exotic one, 
occurs a t a co nstant pressure, and is accompanied by a density ju mp at t he ph ase interface(Baym 



& Chin (119761 )). However, as shown by Glendenning (Glendenning (|l992l ). (ll997l )). the properties of 
phase transition with more than one conserved charge are quite different from the constant pressure 
transition. In such a case, a mixed phase(MP) can be made in the interior NSs. The NSs containing 
MP matter is so-called hybrid stars (HSs) by many investigators. 

It is known that a NS will spin down due to braking (e.g. electric-magnetic radiation or 
gravitational wave radiation). Deconfinement transition can proceed to occur in HSs during the 
spin down. Such evolutionary processes induce not only the changes in stellar structure but also 
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continuous release of latent heat. In the present paper, we calculate the energy release due to a 
phase transition in a rotating HS. The heat arises only at the phase interface for a constant pressure 
transition. In our case, the heat release will be distributed over a MP region. Our calculations are 
based on the perturbation theory developed by IIartle(IIartle (|l967l )). Regard EOS for hadronic 
matter and strange quark matter as the input, we will self-consistently acquire the deconfinement 
heat luminosity with respect to the characteristic age of the star. 

In this work, we choose the simplest possible nuclear matter composition, namely neutrons, 
protons, electrons, and muons (npe/x matter) and ignore superfluidity and superconductivity. 

The paper is organized as follows. In Sec. 2, we introduce notation and describe general prop- 
erties of deconfinement phase transition in stellar with particular emphasis on the existence of the 
MP and varying pressure in the MP. The rotating structure of star with deconfinement phase tran- 
sition are presented in Sec. 3. In Sec. 4 we present the calculation of deconfinement phase transition 
heat luminosity associated with EOS of the MP and structure of rotating star. The conclusion and 
discussions are summarized in Sec. 5. 



2. Deconfinement phase transition 



Quark deconfinement phase transition is expected to occur in neutron matter at densities above 
the nuclear saturation density ps = 0. 16/m ~^. Early works on the possible occurrence of quark 
matter in neutron stars(Baym & Chin (jl976l )) were based on the assumption that hadron matter 
and quark matter were both charge neutral(with only one independent chemical potential). As a 
consequence, the transition was described using Maxwell construction and the resulting picture of 
the star consisted of quark matter core surrounded by a mantle of hadron matter, the two phases 
being separated by a sharp interface. 



In the 90s, Glendenning(Glendenning (jl992l ). (jl997l )) pointed out that this assumption is too 
restrictive. More generally, the transition can through the formation of a MP of hadron matter 
and quark matter, total charge neutrality being achieved by a positively charged amount of hadron 
matter and a negatively charged amount of quark matter. Therefore, at present, most of the 
approaches to deconfinement matter in NS matter use a standard two-phase description of EOS 
where the hadron phase(IIP) and the quark phase(QP) are modelled separately and resulting EOS 
of the MP is obtained by imposing Gibbs conditions for phase equilibrium with the constraint that 
bary on nu mber as well as electric charge of the system are conserved (Glendenning (|l997l ). schertler 
et al (|2OO0l 'l). 

We have to deal with two independent chemical potentials (pmPe) if we impose the condi- 
tion of weak equilibrium. The Gibbs condition for mechanical and chemical equilibrium at zero 
temperature between the HP and the QP reads 



PHp{^J'n,^J'e) = PQpi^J'UjfJ'e 



(1) 



where pnp is the pressure of HP and p q p is the pressure of QP. We use the EOS of the relativistic 
mean field model (Glendenning ()1997l )) for hadron matter and employ an effective mass bag- model 
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EOS for quark matter (Schertle et al. (|l997l )). Only two independent chemical potentials remain 
according to the corresponding two conserved charges of the /3-equilibrium system. The volume 
fraction occupied by QP for every point on the MP curve 

X = ^'^^ (2) 
Vhp + Vqp 

can be obtained by imposing the condition of global charge neutrality in the MP 

Xqqp + (1 - X)<1HP = 0. (3) 

Where q denotes the charge density. Finally, the total energy density e and baryon number ps can 
be calculated using 

e = X^QP + (1 - x)eHP (4) 
PB = XPQP + - X)PHP- (5) 

Taking the charge neutral EOS of the HP, Eq.(l), (2) and (3) for MP and the charge neutral EOS 
of the QP, we can construct the full hybrid star EOS. The resulting of chemical potentials are 
shown in Fig.l. Below the charge neutral HP curve and above the charge neutral QP curve the 
HP is positively charged {qmp > 0) and the QP is negatively charged((7HP < 0). Therefore, the 
charge of hadronic matter can be neutralized in the MP by an appropriate amount of quark matter. 
In the MP the volume proportion of quark phase is monotonically increasing from fin^ to fin^. In 
Fig. 2 we show the model EOS with deconfinement transition which is the typical scheme of a first 
order transition at finite density with MP. The phase transition construction in a two-component 
system leads to continuously increasing pressure of the MP with increasin g deii sity. We choose the 



parameters for hadronic matter EOS which have given by Glendenning (|l997l ) and quark matter 



EOS with s quark mass = 150MeV, bag constant B^^^ = 160MeV, coupling constant g = 3. 



3. Rotating evolution of hybrid stars 



With the evaluated hybrid star EOS presented above we now turn to analy se the structure of 
the co rresponding rotating HSs. Using the Hartle's perturbation theory (|l967l ). Chubarian et al 
(|2000|) have studied the change of the internal structure of the HSs due to rotation. In this paper, 
we also apply Hartle's approach to investigate the structure of rotating HSs. Hartle's formalism is 
based on treating a rotating star as a perturbation on a non-rotating star, expanding the metric of 
an axially symmetric rotating star in even powers of the angular velocity 17. The metric of a slowly 
rotating star to second order in the angular velocity fi, can be written as 



e^^*-) [1 + 2{ho + h2P2)]dt' + e^^^> [1 + 



A(r)r 



2(mo -I- 7722^2) 



]dr'^ 



+r^[l + 2iv2- h2)P2]{d9^ + sin^ 



(r - 2M(r)) 
w{r,9)dtf} + 0{n^ 



(6) 



Here e^^*"), e^^'') and M(r) are functions of r and describe the non-rotating star solution of the 
Tolman-Oppenheimer-Volkov (TOV) equations(Oppenheimer& Volkoff (119391 )). P2 = P2{9) is the 
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l = 2 Legendre polynomials, cu is the angular velocity of the local inertial frame and is propor- 
tional to the star's angular velocity whereas the perturbation functions hQ,h2,mQ,m2,V2 are 
proportional to $7^. we assume that matter in the star is described by a perfect fluid with energy 
momentum tensor 

T"'' = (e + P)u''u'' + Pg"" (7) 

The energy density and pressure of the fluid are affected by the rotation because the rotation 
deforms the star. In the interior of the star at given (r, in a reference frame that is momentarily 
moving with the fluid, the pressure and energy density variation is respectively 

5Pir,e) = [e{r) + P{r)M +P2^'2W] (8) 

6e{r,e) = ^[e{r) + P{r)U+P*2P2{e)] (9) 

here, and dimensionless functions of r, proportional to O^, which describe the pressure 

perturbation. The rotational perturbations of the star's structure are described by the functions 

ho, rn(),pQ,h2,m2,V2,P2- These functions are calculated from Einstein's field equations. The effect 
of rotation described by the metric on the shape of the star can be divided into contributions: A 
spherical expansion which changes the radius of the star, and is described by the functions ho and 
mo. The other part is a quadrupole deformation, described by functions h2, V2 and m2. As a 
consequence of these contributions, the difference between the gravitational mass of the rotating 
star and the non-rotating star with the same central pressure is 

6Mgrav = mo{R) + (10) 
The change in the radius of the star is given by 

5R = UR) + UR)P2{0) (11) 

We wish to study sequences of the rotating stars with constant total baryon number at variable 
spin frequency u = Q/2Tr. The expansion of total baryon numbers in powers of O is 

Nb = N^ + SNb + 0{n^) (12) 

where 

N%= I nB{r)[l-2M{r)/r]-^>'^ATTr'^dr (13) 
Jo 

is the total baryons number of non-rotating star and 

+^!=^i^(. + FW(0}4>r.=* (14) 

here uin is the rest mass per baryon. To construct constant baryon number sequences, we first solve 
the TOV equations to find the non-rotating configuration for giving a central pressure P(r=0). And 
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then, for an assigned value of the angular velocity Q the equations of star structure are solved to 
order Q^, imposing that the correction to the pressure Po(r = 0) being not equal to zero. The value 
of pQ{r = 0) is then changed until the same baryon number as that non-rotating star is obtained. 

The results for the stability of rotating HSs configurations with possible deconfinement phase 
transition according to the EOS described above are shown in Fig.3, where the total gravitation 
mass is given as functions of the equatorial radius and the central baryon number density for 
static stars as well as for stars rotating with the maximum rotation frequency i^fe. The dotted 
lines connect configuration with the same total baryon number and it becomes apparent that the 
rotating configurations are less compact than the static ones. In order to explore the increase 
in central density due to spin down, we create sequences of HSs models. Model in a particular 
sequence have the same constant baryon number, increasing central density and decreasing angular 
velocity. Fig.4 displays the central density of rotating HSs with different gravitational mass at zero 
spin, as a function of its rotational frequency. In the interior of these stars, the matter can be 
gradually converted from the relatively incompressible nuclear matter phase to more compressible 
quark matter phase. 



4. Deconfinement heating rate 

As the star spins-down, the centrifugal force decreases continuously, increasing its internal 
density. Fig.4 identifies the fact that quarks are accTimulating in the interior of the star with 
decreasing rotation frequency u. Since the deconfined phase transition is first order phase transition, 
there are latent heat produced with the transformation of hadron matter into quark matter. The 
deconfinement phase transition heating play an important role in the process of compact star's 
thermal evolution. 

Corresponding to the EOS depicted in Fig.2, Fig.5 shows that energy per baryon of HS matter 
as a function of the baryon number density. Two intersectant solid lines denote the pure hadronic 
matter phase and pure quark matter phase respectively. The deconfinement transition just occurs 
at point 1 until point 2 after which a pure quark matter phase appears. Comparing HP curve to 
MP one, we find that the enthalpy increase in MP is slower than HP when baryon number density 
increases. It is a vivid representation of latent heat release when a phase transition occurs. 

How to describe the latent heat is now a key issue. It is quite evident matter that the energy 
release due to the phase transition is direct proportion to the difference between the HP derivative 
and the MP one at point 1. Similarly, we can express the energy release per baryon as 

5e-6e = {^-^)dpB (15) 

OPB OPB 

where denotes the enthalpy change per baryon for any density in MP curve, density increase 
assumed no phase transition proceeds to occur. For example, point A in Fig.5 has two possible 
enthalpy change ways, AC and AD, when the density increases. No transition occurs along AC 
while AD corresponds to the real situation. In order to derive We rewrite the expression of 
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HP and QP volume Vqp = ^^^,Vhp = and define parameter r/ = Nq/Nb- Substitute these 
equation into Eq.(2), we can get 

X = (16) 

rjpHP + (1 - ri)pQP 

■nPHP + (1 - mpQP 

Replacing Eq.(16) and Eq.(17) in Eq.(4) and Eq.(5), we have energy density for a given r], 

VPHP ^ (1 - V)PQP 



VPHP + (1 - ri)PQP VPHP + (1 - V)PQP 

~ ^ PQPPHP 

VPHP + (1 - V)PQP 
In such a case, We can obtain the energy per baryon as 



(19) 



Furthermore we get 



-z = V^QP + (1 - v)eHP (20) 



de deqp -.denp 

-V^^ + i^-V)^ (21) 



dpB dps dps 
Energy release per baryon at point A in the process of deconfinement phase transition can be given 
using above Eq.(15) and (21). 

Ae = Se - 6e = (,^ + (1 - .)^ - ^)6pB (22) 
dps opB ops 

In Eq. (22), we see Ae equals when 77 takes 1 or 0, returning to pure phases, HP or QP. 
Integrating for whole star, we get the total latent heat release unit time for a star 



H 



de f de „^ , , 



Since the rotation frequency and its derivative can be observed, we can rewrite the Eq. (23) as 

r dp 

H= —0{t)pBdV (24) 

The heat release per baryon || due to reduction in rotating frequency is simulated for a given 
rotating sequence N = I.QQNq (the static mass M = l.SM©), is plotted in Fig.6. We find that the 
energy release effectively is enhanced with increasing density and frequency. So the deconfinement 
transition in the star is stronger during the early ages, and the energy release rate becomes larger 
and larger from the outer region to the center of the star. 

Of course, the heat luminosity H can be obtained only if the spin-down rate of a star given. 
The common case is one induced by magnetic dipole radiation. It reads 

^ = _i^^Vsin2 (25) 
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is induced by magnetic dipole radiation, where I is the stellar moment of inertia, /i = \BI^ is the 
magnetic dipole moment, and 9 is the inclination angle between magnetic and rotational axes. 

The total heat luminosity as a function of the characteristic age for different magnetic fields, 
for the rotating sequence = l.GGA^© (the static mass M = I.SMq), is presented in Fig. 7. 
Heating source inside the stars is an important factor effect on the cooling of NSs. NSs first 
cools (for t < lO^yrs) via various neutrino emission before the surface photon radia tion t a ke ove r. 
Several heating mechanisms, for example, rotoc hemica l heating(Reisenegger et al. (119951'). (l2006l )). 
compositional transitions in crus t (lid a & Sato (|l997l )). crust cracking (Cheng et al. (jl992l )) and 
vortex pinning (Van Riper et al. (|l995l )). have been discussed in detail. It is generaly expected that 
the heating sources significantly contribute to surface emission for the old NSs with low magnetic 
fields, characteristic of millisecond pulsars. Recently the thermal emission data showed to have 
high surface temperature for a few of millisecond pulsars. Especially for PSR J0437-4715, the 
temperature can be inferred as 1.2 x IQ^K. In light of model f or NS thermal evolution including 
rotochemical heating source, Reisenegger (Reisenegger et al. ([200a)) found that the theoretical 
result is 20% lower than the inferred temperature. It is noteworthy in Fig. 7 that our heat luminosity 
for magnetic field ~ lO^G and ~ lO^G lasts during a term far longer than IQ^yrs and is much higher 
than the other heat generation. We thus think that high temperature of some millisecond pulsars 
with low magnetic fields (Kargaltsev et al. (|2004l )) can be explained using predicted heating 
model of HSs. For PSR J0437-4715, we predict the surface emission L^oi = Hiv,!))^]. — ^^J^ ) 
with rotation period and its derivative P = 5.74:, P = 3.64 x 10^^^. The result is estimated 
as Lboi ~ 6.6 X 10^^ergs~^ considerably consistent with the observed thermal X-ray luminosity 
Lx = iiraR'^T^ ~ 2.5 x W'^^ergs'^. 



5. Conclusions and discussions 

The nuclear matter can continuously be deconfined to quark matter during the spins-down of 
star. The deconfinement phase transition heating of rotating hybrid stars have been investigated in 
this work. The latent heat release of such a first order phase transition associated with the change 
of energy per baryon during transition and rotational structure evolution of star. Using Hartle's 
perturbative approach, we have calculated the change of internal structure of rotating hybrid stars. 
For a process of infinitesimal compressing during the rotation of star, we can get the energy release 
per baryon by inputing the EOS of the MP. Furthermore, we self-consistently get the total energy 
release combining with the rotational structure of star. The results show that the latent heat release 
per baryon per unit frequency enhances with the increase of baryon number density and rotational 
frequency in HSs. 

We have calculated the changes in the internal structure of a compact star during its spin-down. 
As shown in previous studies the deconfinement phase transition occurs if quark matter is stable 
state at high densities, since the nuclear matter is compressed in the interior of the star during 
spin-down. We consider a first-order transition presented by Glendenning (1992) to calculate the 
release of latent heat when the star is slow down. The heating rate arising from the processes has 
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been estimated. We find its significant efi'ects on the thermal evolution and the effects is much more 
important than the past present heating mechanisms. For old neutron stars with low magnetic fields 
of ~ 10® — lO^G, especially, our predicted temperature is nearly identified with the values observed 
from millisecond pulsars. In future, we expect more observational examples in investigating effects 
of deconfinement heating mechanism on the NSs thermal evolution. Thermal evolution curves need 
to be compared to more data. In this paper, for si mplicity we neglect Coulomb and surface effects 



on the MP. Many investigations (Endo et al (120061 ) )have identified that both effects would restrict 



the region of the MP in the core of hybrid stars, which can effect the release of latent heat. Another 
problem which remains to be investigated is the unified description of middle-age and old pulsars. 
HS model may be no bad selection when our combining deconfinement heating, superfluidity effects 
in nuclear matter together. 

This work is supported by NFSC under Grant Nos. 10603002 and 10373007. 
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Fig. 1. — Electron chemical potential /ig as a function of the neutron chemical potential The 
HP EOS is a rclativistic mean-field model, the quark matter is effective mass MIT bag model with 
rus = IbOMeV, B^/^ = leOMeF, coupling constant g = 3.0. 




Fig. 2. — Model EOS for the pressure of hybrid star matter as a function of the baryon number 
density. The HP EOS is a relativistic mean-field model, the quark matter is effective mass MIT 
bag model with = IbQMeV , B^/'^ = 160MeF,coupling constant g = 3.0. 




Fig. 3. — Gravitational mass M as a function of the equatorial radius (left figure) and the central 
density (right figure) for rotating hybrid stars configurations with a deconfinement phase phase 
transition. The solid curves correspond to static configurations, the dashed ones to those with 
maximum rotation frequency i//.. The lines between both extremal cases connect configurations 
with the same total baryon number. 




Fig. 4. — Central density as a function of rotational frequency for rotating hybrid stars of different 
gravitational mass at zero spin. All sequences are with constant total baryon number. Dash 
horizontal lines indicate the density where quark matter is produced. 
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Fig. 5. — The energy of per baryon as a function of the baryon number density for hybrid star 
matter. Two intersectant sohd hnes denote the pure hadronic matter phase and pure quark matter 
phase respectively. Line AC represents the state of no deconfined phase transition during rotation 
of star. 
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Fig. 6. — For 1.5 Mq hybrid star, the energy release of per baryon unit frequency as a function of 
the baryon number density in mixed phase region for various frequency. 
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Fig. 7. — For 1.5 Mq hybrid star, the dcconfinemcnt phase transition heating rate change with 
rotational evolution of star for various magnetic fields. 



